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Page 4: The second and third paragraphs of the APPLICATIONS section should 

be replaced by the following paragraphs: 

Reference 13 suggests that for turbulent boundary- layer separation 
H^ ranges from 1.8 to 2.6. With E±^2 = 2.2^ the Mach number ratio M 2 /M 1 
obtained from equation (7) is 0.762. As shown in figure 3(a)^ this Mach 
number ratio yields pressure ratios from the oblique-shock relations that 
are in good agreement with ’’first peak” pressure ratios obtained on forward- 
facing steps in references 8 to 11. The data of reference 12 fall below the 
rest; however, figure 9 of this reference suggests that these data may have 
been taken in the transition region. A Mach number ratio of 0.81 yields 
pressure ratios in good agreement with ’’inflection-point” pressure ratios 
obtained on wedges (fig. 8 of ref. 8). According to reference 8, these 
first-peak- and inflection-point pressure ratios provide first approximations 
to the pressure ratios for which separation is likely to be encountered on 
control-surface configurations. 

The separation pressure ratios obtained from the obliq ue-sho ck relations 
and equation (4) of reference 4 for a Mach number ratio of VO.494 = 0.703 
are presented in figure 3(b) as are also the pressure ratios obtained from 
equation (8) of reference 4. These pressure ratios are compared, as in 
reference 4, with the ’’inflection point” pressure ratios for wedges reported 
in reference 8. Inasmuch as the results from equations (4) and (8) of 
reference 4 should be in agreement, the discrepancy shown in figure 3(b) 
suggests that the oblique-shock linearization used in equation (8) of 
reference 4 is inadequate. 

Page 12: The sentence following equation (Cl) should read ’’Using relations 

(B4)” instead of ’’Using relations (B2)”. 

Page 20: Figure 3 should be replaced by the attached figure. 
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(a) First peak pressure ratio for 
forward-facing steps. 


(b) Inflection-point pressure ratios 
obtained on wedges . 


Figure 3. - Pressure ratios associated with shock-induced, turbulent, boundary- 
layer separation. 
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TECHNICAL NOTE 3454 

EFFECT OF A DISCONTINUITY ON TURBULENT BOUNDARY-LAYER -THICKNESS 
PARAMETERS WITH APPLICATION TO SHOCK-INDUCED SEPARATION 
By Eli Reshotko and Maurice Tucker 


SUMMARY 

The problem of shock-induced tirrhiilent boundary-layer separation 
was analyzed by an approximate method. This method is based on a moment - 
of -momentum eqmtion and calculates the change of boundary-layer-thickness 
parameters and form factor caused by a discontinuity where the effects of 
friction can be neglected. The form of the result suggests that the Mach 
number ratio across the shock is a characteristic parameter for defining 
shock-induced separation. The method is also lised to estimate the effects 
of mass transfer on the boundary-layer-thickness parameters for zero pres- 
sure gradient. 


INTRODUCTION 

Because shock- induced separation of the boundary layer is of partic- 
ular interest^ a nxmxber of recent analytical and experimental studies have 
been concerned with this phenomenon. Owing to the steep pressure grad- 
ients involved when a shock wave interacts with a turbulent boundary layer, 
the most promising analytic attacks on the interaction problem have been 
"discontinuity" analyses such as those of references 1 to 4 wherein effects 
of friction are considered negligible compared with effects of pressure 
gradient. Tyler and Shapiro (ref. 2) have ass\maed a power-law velocity 
profile ahead of the shock and a "separating" velocity profile downstream 
of the shock with constant pressxrre gradient between these stations. 

Crocco and Probstein (ref. 3) have applied the mixing theory of reference 
5 in their analysis. Mager (ref. 4) has utilized the semi-en^Jirical 
Gruschwitz auxiliary equation in conjunction with the Stewartson transfor- 
mation (ref. 6) to predict the pressure rise required for shock- induced 
turb\ilent boundary -layer separation. 

The present analysis is concerned with the effect of a discontinuity 
on a nonuniform velocity -profile characteristic of a turbulent boundary 
layer. In this analysis, a discontinuity is considered as an abrupt but 
not necessarily discontinuous change in the flow field. Although friction 
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is not to be considered^ it is convenient to describe the discontinuity- 
effects in terms of conventional bomdary-layer-thickness parameters . 

The analysis utilizes a moment -of -momentum equation similar to that of 
reference 1 , but, of co\rrse, without the shear terras. The rational 
nature of this moment -of -momentum equation, which retains the same form 
for both two-dimensional and axially symme-trlc flows, permits its appli- 
cation to a number of flow problems. For shock- induced separation, re- 
sults similar to those of reference 4 are obtained and are compared with 
the experimental results of references 8 to 12. The effect of a discon- 
tinuity involving mass transfer (blowing or suction) upon the thickness 
parameters is also indicated for the case of zero pressure gradient. 


ANALYSIS 


The techniques of turbulent bo\mdary-layer analysis are utilized for 
the present analysis. With the assumption of a constant stagnation tem- 
perature transverse to the flow siirface, consideration of the rate of 
change of momentum in the longitudinal direction leads to an integral 
equation identical -with the Karman momentum Integral equation except 
for absence of the shear terms. The corresponding moment -of -momentvim 
equation is then obtained through multiplication of the integrand of the 
momentum Integral equation by a distance normal to the surface Y and 
then integrating with respect to Y. (All symbols are defined in appendix 
A; appendix B outlines the development of the moment -of -momentum eqimtion 
with provision for mass transfer using a form of the Stewartson transform- 
ation. ) 


The transformed moment -of -momentum equation with mass transfer absent 
takes the form 

dUe H.(h 2 - l)(H. + l) 


<^1= - U 


( 1 ) 


Use of the transformed equation simplifies the analysis in that effects 
of Mach number need not be considered in treating the thickness param- 
eters. Thus, takes on values typical of the form-factor values en- 
countered in incompressible -flow treatments. When the transformation 
relation = Ma^ is used, equation (l) may be written 


M 


2dH^ 

H.(h2 - 1)(H. + 1) 


( 2 ) 


Upon integration, equation (2) becomes 


AM = 


H? l/(Hi+l) 

(H^ - l)^^^(% + 1) 


= f(Hi) 


( 3 ) 
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where A is an integration constant. The right-hand side of equation 
( 3 )^ which will he denoted as f(Hj^)^ is plotted in figure 1. Shown for 
comparison in figure 1 is the equivalent of f(Hj^), as proposed in refer- 
ence 4^ which is where is the Gruschwitz form factor. The 

two functions differ only slightly. 


For two-dimensional flow, the change in momentum thickness caused by 
a discontinuity is obtained by using the momentijm integral equation from 
appendix B, again neglecting the mass-transfer term: 


d0n 


dU^ 


(2 + Hi) ^ = 0 


(4) 


Using the relation U^ = MaQ and equation ( 2 ) changes equation (4) to 


d 0 i 2(2 + H^)dHi 

®i " H.(H? - 1)(H. + 1) 


which upon integration yields 


(hJ - 1) ^ (Hi + 1) -l/(Hi+l) 

B0i = -i 5 e 

Hi 


= g(Hi) 


( 6 ) 


where B is an integration constant. The right-hand side of equation 
(6) denoted as g(Hi) is plotted in figure 2. 

Thus, across a discontinuity, the Mach number ratio and momentuun- 
thickness ratio from eqi,iations (3) and (6) are, respectively 


Mg 



( 7 ) 

( 8 ) 


where the subscript 1 refers to conditions aiiead of the discontinuity 
and the subscript 2 refers to conditions behind the discontinuity. 
Eq\iation (?) gives the magnitude of the discontinuity in terms of Mach 
number ratio necessary to change the transformed form factor from 

to H. o and can also be thought of as the relation that gives the change 

c, 

in form, factor resulting from a discontinuity of arbitrary but known mag- 
nitude. In applying equations (?) and (s) and figures 1 and 2 to com- 
pressible flow^ the following relations between compressible quantities 
and their transformed equivalents may be useful: 
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+ 




0 




_r+i 

,2^2(r-i} 


( 9 ) 

( 10 ) 


These relations can he obtained from the conventional definitions of the 
boundary-layer thickness parameters H and 0 as shown in appendix C . 
The effect of a discontinuity on displacement thickness follows from the 

definition 


APPLICATIONS 
Shock-Induced Separation 

Equations (?) and (8) may be applied to the problem of shock-induced 
separation of the turb\lLent boundary layer if appropriate values of 
are selected. For zero pressure gradient ahead of the shock wave a form 
factor H. = 1.286 is assumed, which corresponds to a transformed 
seventh-power- law velocity profile. Such an assumption yields form fac- 
tors for the physical compressible flows that are compatible with exper- 
imental observations. However, as indicated by equation (Cll) of appen- 
dix C, the physical profile is not a power-law velocity profile. 


Reference 13 suggests that for turbulent bo\andary-layer separation 
ranges from 1.8 to 2.6. With = ^.2, the Mach number ratio 


Hi „ _ 

obtained from equation (?) is 0.762. As shown in figure 3, this 

Mach number ratio yields press\ire ratios from the oblique-shock relations 
that are in good agreement with "first peak" pressure ratios obtained on 
forward-facing steps in references 8 to 11. The data of reference 12 
fall below the restj however, figure 9 of this reference suggests that 
these data may have been taken in the transition region. A Mach number 
ratio of 0.81 (obtained with H^^g = yields pressure ratios in good 

agreement with "inflection-point" pressure ratios obtained on wedges (fig 
8 of ref. 8). According to reference 8, these first-peak- and inflection 
point pressure ratios provide first approximations to the pressure ratios 
for which separation is likely to be encountered on control- surface 
configurations . 

The pressure ratios obtained from the oblique-shock relations and 
equation (4) of reference 4 which states, in the notation of this paper, 

that — = -,70.494 = 0.703 are presented in figiire 3 as are also the 
M]_ ^ 

pressure ratios obtained from equation (8) of reference 4. Inasmuch as 
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the results from eqmtions (4) and (8) of reference 4 shoiild be in agree- 
ment^ the discrepancy shown in figure 3 suggests that the oblique-shock 
linearization used in equation (8) of reference 4 is inadequate. The 
resiilt of reference 4 overestimates the pressure ratios obtained experi- 
mentally on forward-facing steps. 

Effects of pressure gradient . - The present analysis can be applied 
to the case of favorable or adverse pressure gradient ahead of the shock 
when the appropriate form factors are known. Qualitatively it is known 
that the fonn factors for adverse and favorable pressiore gradients are, 
respectively, greater and less than the flat -plate form factor. Thus, 
for favorable pressure gradients ahead of the shock, a stronger shock is 
required, whereas for adverse pressirre gradients ahead of the shock, a 
weaker shock would separate the turbulent boundary layer. 

For the case of boundary-layer development ixnder an adverse pressure 
gradient large relative to the skin-friction term but where the incident 
shock is of insufficient strength to separate the flow, equations (7) to 
(lO) may be used to obtain the change in H^, 0^^, and 5j^ across the 
shock. The remaining development to separation requires solution of the 
equations of the compressible turbulent boimdary layer with a stipilated 
skin-friction relation. 


Effect of Mass Transfer 


The effect of suction or blowing on form factor for two-dimensional 
flat -plate flow may be approximately evaluated by using the momenta 
integral and moment -of -momentum equations of appendix B. The method 
requires that the skin-friction coefficient be negligible relative to 

the mass-transfer coefficient In the vicinity of the mass -transfer 

f’e’^e 

region, the press\ire gradient terms may be appreciable. The following 
analysis, which omits consideration of pressxure gradient, may thus yield 
only qualitative results. A similar analysis for flow over a cone is 
given in appendix D. For two-dimensional flow, the moment-um integral and 
moment -of -momentum eqmtions are, respectively. 


and 


f£i 

dx ■ Ug 


0 




1 ) — 

Ue 


( 11 ) 


( 12 ) 


Eliminating the quantity 
yields 




^i 

H? - 1 


between equations (ll) and (l2) 



0 . 

1 


6 


MCA TN 3454 


which upon integration becomes 

ce. = 


- 1 


Hi + 1 




where C is an integration constant, or 


0 

0 


i^ 

i.l 



(13) 


(14) 


where the subscripts 1 and 2 refer, respectively, to conditions be- 
fore and after the mass transfer. The function R(Hj^) is plotted in 
figure 4. 

The relation between the transformed momentum-thickness ratio and 
the mass transfer is obtained from the momentum equation with the assimip- 
tion that V^/Ug is constant for the region of mass transfer. Thus upon 

integration, equation (ll) becomes 


0 .. 


PwV 


1.2 p^a^M ®i,l 


or, since m = represents the mass flow per \anit width 


= 1 + 

®i,l PQ^O^i,! 


(15) 


Expressing the mass transfer as a fraction of the boundary-layer 
mass flow ahead of the transfer region transforms equation (l5) into 


0 

0 


hi 

lA 


1 + 



where the mass-transfer coefficient p is defined as 


P = 


m 

Pe^e^^l 



The quantities % 

law profiles in reference 14. 


(16) 


(17) 


are tabulated for certain typical power- 
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The amount of mass removal by suction ahead of a normal shock that 
is required to prevent flow separation behind the shock may be estimated 
in the absence of pressure gradient ahead of the shock and for negligible 
skin friction. The form factor ahead of the normal shock for incipient 
separation is found from equation (?) and the mass removal required to 
obtain this form factor is determined from equations (14) and (16). 

For an assumed seventh-power-law transformed- velocity profile and with 
Hi 2 = the percentage of boundary -layer mass flow to be removed 

rises steeply with Mach number^ reaches a peeik of about 8 percent at 
M » 3, and gradually decreases at higher Mach numbers. The mass flow 
to be removed, however, will increase with Mach number for constant 
ambient pressure and temperature. In the absence of exp)erimental data, 
these results must be considered qualitative. The method cannot be 
used to estimate the effects of mass transfer behind shock-induced 
separation, since the analysis does not apply after the onset of 
separation. 


A method has been devised for calculating the change of boundary- 
layer- thickness parameters and form factor caused by a discontinuity in 
the absence of friction effects. The form of the result suggests that 
the Mach number ratio across the shock is a characteristic parameter for 
defining shock- induced separation. The experimental data for forward- 
facing steps and for wedges are well described by the curves of constant 
Mach number ratio (ratio of Mach number ahead of to that behind discon- 


also been applied to estimate the amount of mass removal that is required 
ahead of a normal shock to prevent flow separation behind the shock. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, April 11, 1955 


SUMMARY OF RESULTS 


tinuity) 



The method has 
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APPENDIX A 
SIMBOLS 


The following symbols are used in this report: 
Af C integration constants 
a velocity of sound 


f(Hi) 


g(Hi) 


H 

M 

m 

N 

P 

r 


R(H.) 

t 

U 

u 

V 

V 

X 


Hf 

(h| - 1)^/^ (H. + 1) 

(nf - (H^ + 1) ^-i/(%+i) 


form factor; H h S*/0 
free -stream Mach number 
mass transfer per unit width 

power-law velocity profile parameter; 



pressTire 

radius of axially symmetric body 



temperatixre 

transformed longitudinal velocity 
longitudinal velocity 
transformed normal velocity 
normal velocity 

transformed longitudinal coordinate 


CG-2 
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X 

Y 

y 

r 

6 




T1 

M* 

G 


longitudinal coordinate 
transformed normal coordinate 
normal coordinate 
ratio of specific heats 
boimdary- layer thickness 
displacement thickness 
transformed displacement thickness^ 



Gruschwitz form factor^ 1 


^ - 1 

+ 1)J 


mass-transfer coefficient 
momentum thickness 
transformed momentum thickness, 




dY 


dY 


p dens ity 

T shear stress 

Subscripts ; 

e free stream ("external") 

i transformed or "incompressible" value 

s stagnation value 

w wall value 

0 free-stream stagnation value ahead of discontinuity 

1 conditions ahead of discontinuity 

2 conditions behind discontinuity 
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APPENDIX B 


DEVELOPMENT OF MOMENT -OF -MOMENTUM EQUATION 


The behavior of nonuniform profiles of the boimdary-layer type in 
a compressible flow can be studied from continuity and momentum consider- 
ations. In examining this behavior^ the effects of shear are neglected, 
although presumably the nonuniform profiles might have been set up 
through shearing action. With the stagnation temperature in the trans- 
verse direction considered to be constant, the equations of motion are 


Continuity: 


Moment m: 


S (P«) + (P'') = 0 ^ 


Su Sp 


> 




(Bl) 


Equations (Bl) resemble the boundary-layer equations except that the 
friction term is lacking from the momentum equation. The application of 
a Stewartson transformation to equations (Bl) and reduction of the result 
to an integral equation follow the procedure of reference 15. A modi- 
fied Stewartson transformation can be described by the following coor- 
dinate transformations : 


Pe ^•e 
dX = — — dx 
Pq ^0 


dY = ^ dy 

Pq ^0 


(B2) 


The transformed coordinates are now represented by upper-case let- 
ters (X, Y) and the subscript e refers to local conditions at the outer 
edge of the boundary-layer-type profile. The subscript 0 denotes refer- 
ence values, which are here taken as free -stream stagnation values ahead 
of the discontinuity. The relation between transformed and physical lon- 
gitudinal velocities resiilting from the Stewartson transformation is 
a© 

U = — u or for the external flow, U^ = Muq. 
e 

From the method of reference 15, the integrand of the momentum in- 
tegral eqmtion is formed. The resvilt for the case of the present paper 
is 


(U(Ue - U)]^ + [V(Ug - U) ]^ + Ue^(Ug - U) = 0 


(B3) 
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When equation (B3) Is integrated with respect to Y between Y = 0 
and Y = where 6^^ is the transformed bovindary-layer thickness, the 
transformed momentim integral equation is obtained: 


d0i 0i(2 + Hi) dUg 

dX Ug dX ■ Ug 


(B4) 


where V^/Ug is the term representing the effect of mass transfer. 

If the integrand (B3) is multiplied by Y and the resxilting rela- 
tion integrated with respect to Y, the moment -of -momentum relation is 
obtained as indicated, for example, in reference 7. The moment-of- 
momentum equation in the transformed plane can be written 


e 


dH. 
1 

i dX 


0. dU H.(H? - 1)(H. + l) 
1 e 1 * 

Ug dX 2 




(B5) 


Application of equations (B4) and (B5) to tm^bulent flow implies 
that densities and velocities are time-averaged quantities and that the 
correlation terras (those involving products of pertirrbation q\iantities) 
are negligible. 
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APPENDIX C 


RELATION OF COMPRESSIBLE AND TRANSFORMED 
BOUNDARY-LAYER FORM FACTORS 

The relation between compressible and transformed boundary-layer 
form factors will be formulated through the individml relations for 
momentum and displacement thicknesses^ respectively. Although the pro- 
cedures of the present paper are for flows with constant stagnation tem- 
perature in the transverse direction, the results of this appendix apply 
also to the case of variable stagnation temperature in the transverse 
direction. 


Momentum Thickness 


The momentum thickness of a compressible boundary layer is defined 


as 


e = 


^0 


-en-fi 

Pe% \ ^e) 


dy 


(Cl) 


Using relations (B4) and the knowledge that = ~ changes equation 
(Cl) to 


Ue ^e 


where 




dY 


(C2) 


(C3) 


Displacement Thickness 


The displacement thickness of a compressible boundary layer is de- 
fined as 

_6 ^ 

' dy (C4) 

, \ / 

/o 


6 * = r (i . 

Jo \ "e^ej 
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With constant pressure through the ho\mdary layer, equation (C4) can be 
rewritten as 

The static -temperatiire distribution in the boundary layer can be written 



(C6) 


where t is the stagnation temperature corresponding to the static 
s 

temperature t . 



which for constant stagnation ten 5 )erature through the boimdary layer is 



Since H = 6*/0 and eqiiations (C2) and (C7) with the 

definitions (C3) and (C8) yield 

H = ^ ~ - + l) (CIO) 


or 




M 


( 6 ) 


14 


MCA TN 3454 


The Stewartson transformation causes distortion of the velocity pro- 
file with increasing Mach number M. Assume that the transformed- and 

/ \Vn 

physical -velocity profiles are the power-law profiles -^ = and 


■© 


1/Ni 


_U 

U, 


leads to the relation for 


With — = 

^e % 


and 

r = 1.4 


‘f 


Pq^O 

pa„ 


dY, the transformation 



5 + 


2 + NiVuJ 


5 + 


2 + 


(Cll) 


Eqmtion (Cll) indicates that if is taken to be constant as in a 

standard power-law velocity profile_, then the physical profile parameter 
N will vary with the normal coordinate y. This variation is accentuated 
with increasing Mach number. 
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APPENDIX D 

DERIVATION OF RELATIONS FOR AXIALLY SYMMETRIC FLOW 

For axially symmetric flow, the equations corresponding to (B4) and 
(B5) are 

Momentum: 


d6i 0i(2 + Hi) dUe 9i 
dX"^ Ug dX‘*'rdX"Ug 


Moment of Momentum: 
dX " 


V, 


(h| - 1) - 


Sj dUg Hi(Ef - 1)(% + 1) 

U„ dX 2 


(Dl) 


(D2) 


The moment -of -momentum equation is the same for axially symmetric 
as for two-dimensional flow. Thus, 


M2 _ f(Hi,2) 


(7) 


is correct for axially symmetric flow. Examination of the momentim inte- 
gral equation (Dl) shows that the axially symmetric equivalent of equa- 
tion (s) is written 


(r0i) g(Hi 2 ) 

H = — (D3) 

(r0i)i g(Hi^i) 

Sitice across a discontinuity, the radii r^ « T 2 , equation (D3) is es- 
sentially that for two-dimensional flow. 


Consideration of suction or blowing on a cone in a manner analogous 
to that for two-dimensional flat-plate flow leads to the following equa- 
tions as equivalents for equations (l4), (l5), and (I 6 ), respectively: 


(r0i) 


2 






1 + 



(D4) 

m 

(D5) 
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- H. 


(D6) 


where m now represents the total mass transfer and which is still 
the ratio of mass transferred to the boundary-layer mass flow^ is now 
defined 
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Figure 1. - Function for determining Mach number ratio. 
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Mach nimber ahead of shock, 


Figure 3. - Pressure ratio for shock- induced, 
turbulent, boundary-layer separation. 


H. +1 


NACA TN 3454 


21 



r 


NACA - Langley Field, Va. 



